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ABSTRACT 


Range  and  angle  errors  produced  by  radio  ray  refraction  are  shown 
to  be  calculable  from  closed  form  exponential  expressions.  These  expressions 
are  developed  from  the  error  function  equations  derived  by  Freeman  and 
apply  to  propagation  in  an  exponential  model  atmosphere.  The  exponential 
expressions  are  further  empirically  modified  to  provide  error  data  using  the 
directly  available  radar  tracking  parameters,  that  is,  apparent  elevation 
angle  and  range.  It  is  further  shown  that  the  elevation  and  doppler  error 
angles  may  be  derived  from  the  exponential  expressions. 
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1. 


INTRODUCTION 


The  correction  for  range  and  angle  errors  during  real-time  radar  or 
optical  tracking  operations  generally  requires  the  piecewise  solution  of 
Snell's  ray-tracing  equations  using  digital  computers  Two  approaches  can 
be  taken. 

In  the  first  case,  the  atmosphere  is  frozen  in  time  and  space  and  the 
various  errors  are  calculated  based  upon  a  description  of  the  environment 
obtained  from  direct  measurements  with  instrumentation.  During  the  actual 
tracking  operation  the  measured  ranges  and  angles  are  then  corrected  by 
extracting  the  error  data  from  storage  This  procedure  may  be  impractical 
due  to  the  assumption  of  a  frozen  atmosphere,  instrumentation  limitations, 
the  availability  of  sufficient  computer  storage  and  the  data  access  time 

In  the  second  case,  the  requirement  to  measure,  assimilate,  and 
process  environmental  data  is  removed  by  the  tacit  assumption  that  the  en¬ 
vironment  can  be  defined  by  a  simple  analytical  model,  for  example,  an  ex¬ 
ponential  decrease  of  the  index  of  refraction  with  height.  This  definition  may 
be  quite  acceptable;  in  which  case,  the  propagation  errors  can  be  more 
readily  determined  in  real-time.  Even  in  this  simplified  situation  the  com¬ 
putation  time  required  to  calculate  the  errors  from  the  ray  tracing  equations 
is  not  insignificant,  if  a  high  calculation  accuracy  is  to  be  maintained. 

It  is  therefore  advantageous  to  obtain  simple,  closed  form  equations 
to  calculate  tracking  errors  in  real-time.  This  approach  eliminates  the 
requirement  for  ray-tracing  calculations,  storage  of  error  data  and  avoids 
limitations  in  the  access  time  to  retrieve  stored  data. 

Since  the  index  of  refraction  tends  to  decrease  exponentially  with 
height,  the  exponential  model  has  gained  wide  application.  Studies  in  its 
use  have  led  to  the  well-known  CRPL  surface  corrected  model  which  attempts 
to  adapt  itself  to  locally  varying  conditions  1,s 

Using  an  exponential  model  Thayer  and  Freeman  showed  that  errors 
produced  by  refraction  could  be  represented  by  error  function  expressions. 

This  paper  shows  how  these  expressions  may  be  replaced  by  ex¬ 
ponential  functions  and  emp  irically  adjusted  for  use  at  low  elevation  angles 
and  all  practical  ranges  Finally,  these  exponential  expressions  are  used 
to  derive  equations  to  calculate  elevation  and  doppler  error  angles 


1 


In  the  development  of  these  expressions  the  tracking  parameters, 
where  introduced,  are  in  terms  of  the  apparent  elevation  angle  and  range  of 
the  target  since  these  variables  are  directly  available  from  tracking  systems. 
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2. 


THE  SIMPLIFICATION  AND  MODIFICATION  OF  THE  EQUATIONS 
DERIVED  FROM  SLANT -PATH  INTEGRATION 


The  solution  of  the  range  error  and  total  path  bending  integrals  was 
carried  out  in  Freeman’s  development  by  integration  along  the  slant-path, 

R0  ,  to  the  target.  4,6  (Reference  Figure  1  .  )  This  approach  is  justified  by  the 
fact  that  the  range  error  due  to  the  decrease  in  propagation  velocity  is  much 
larger  than  the  increased  geometrical  range  along  the  curved  path  R  It  is 
apparent  that  the  height  to  any  point,  P,  on  path  Ro  can  be  related  to  the 
geometric  distance,  S,  and  the  true  elevation  angle,  /?(, 

Using  an  exponential  refractivity  model, 

N  =  Ns  exp  (-ch)  (1) 

where  Ns  =  the  station  refractivity  (N  units) 
c  =  the  decay  constant  (km  L) 


Freeman  then  showed  that  the  range  error  and  total  path  bending  in¬ 
tegrals  could  be  solved  to  give 


AR  =  10'8  r  N  ds  =  -‘''-■A  eXp.  (gJ  Ng  [erf  (R 
J  10c  sin  Ph 


:0  cos  %  (A—)  +  g) 

^  rn 


o 

erf  (g)] 


(2) 


T  =  -  1 0"e  JNl"  cos  0  dN  =  . 

$0  f  ^0 


1  u  tan  Pc 


[erf  (R0  cos  0O  (A~ )L  g)  -  erf  (g)] 

i.  r  0 


(3) 


where  g  =  (--* 9  )  "  tan  P0 

r0=  the  earth-xadius  (km) 


Again,  the  range  error  used  throughout  this  discussion  is  defined  to 
be  the  difference  between  the  radio  path  length  along  R  (propagation  dependent) 
and  the  geometric  slant  range,  R0. 


The  first  step  in  the  modification  of  these  equations  is  to  note6  that 


1 _ 

X  +  (x2  +  2)  ¥ 


oo 

^  exp  (x2)  J  exp  (-ts) 
x 


dt  < 


1 


X  +  (x2  +  4/TT)J 


(4) 
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FIGURE  1  RAY  GEOMETRY 
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Under  the  condition  that 


>>2 


oo 


then  exp  (x  )  J  exp  (-t2)  dt  ~  1/2  x 


(5) 

(6) 


x 


and  since 


oo 


srf  x  1  -  - exp  (-t2)  dt 

^  X 


(7) 


substituting  for  the  error  function  in  Equations  (2)  and  (3)  and  making  use 
of  the  Equations  (6)  and  (7)  gives 


AR 


N« 


10c  sin  j30 

No 


[1  - 


g  exP  [g2J 


(k  +  g)  exp  [(k  +  g)  ] 


] 


~  10"  tan  f$0 

where,  from  (2) 


[i  - 


g  exp  [gg] 


(k  +  g)  exp  [(k  +  g)2] 


(8) 


(9) 


k  =  R  (t— - — )~cos  P0 
2  r  o 


and  g  =  -  )^tan  &0>>(2)' 


c  _ a  ,  ,c  r0^_  a 


and  (k  +  g)  =  R  (- - )  ‘  cos  £0  +  (— t-2-  )  tan  Pn  >>(2) 

2  r0  2 


(10) 

(ID 

(12) 


As  the  range  becomes  very  large,  k  increases  and  Equations  (8)  and  (9) 
reduce  to 


Ar 

-  .  NP  _ 

1 0^  c  sin  j30 

(13) 

Ns 

T  = 

1  0e  tan  8n 

(14) 

For  Qo>200  mr,  then  $0  may  be  replaced  by  0O  without  serious  error,  and 
Equations  (13)  and  (14)  then  take  on  a  familiar  form  for  high  angle  tracking.  7 
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A  CORRECTION  AT  SMALL  ELEVATION  ANGLES 


It  was  noted6  that  the  original  Equations  (Z)  and  (3)  give  larger  values 
of  the  range  error  and  path  bending  than  is  actually  true  because  the  slant 
path,  Rq,  lies  below  the  ray  path,  R  Since  the  magnitude  of  refractivity  is 
larger  along  R0  than  at  the  higher  level  along  R,  this  result  is  not  unexpected. 
Replacing  $0  by  the  apparent  elevation  angle,  90,  generates  a  path  which  will 
lie  above  the  true  propagation  path  therefore  giving  values  which  are  too 
small  These  differences  become  most  noticeable,  of  course,  at  small  ele¬ 
vation  angle  s . 

As  the  true  elevation  angle,  $0,  becomes  smaller  the  conditions  on 
the  substitution  for  the  error  function  becomes  invalidated  (reference  Equa¬ 
tion  (11))  In  fact,  it  is  found  that  even  if  is  r  eplaced  by  0O,  the  calculated 
values  from  Equations  (8)  and  (9)  will  always  be  too  large,  particularly  at 
small  elevation  angles 

In  older  to  maintain  the  usefulness  of  the  exponential  functions  at 
small  angles  a  substitution  for  sm  jS0  was  selected  to  be 

sin  j30  =  sin  y]  =  sin  0O  +  k0  (exp  [-kx  Q02])  (15) 

The  constant,  klf  is  selected  to  make  j30  appi  oac  h  0O  at  higher  elevation 
angles  which  has  physical  meaning,  Secondly,  in  the  region  of  zero  elevation 
angles  the  Equations  (8)  and  (9)  remain  finite  and  can  be  constrained  to  pro¬ 
vide  meaningful  results  by  the  suitable  selection  of  k0.  Comparisons  with 
ray-tracing  data,  as  will  be  shown,  indicate  this  form  is  acceptable  and  is 
convenient  to  use  in  calculations 


Equation  (8),  for  example,  then  becomes 


AR 


Nc 


r,  g  exP  [g3] 

L  / 1 ^  i  \  _  Tii- 


i.0  c  sin  X’  (k  +  g)  exp  [(k  +  g)c] 


where  g  and  k  are  defined  as 


g  -  (  2“ ) 2  tan  V 


k  =  R  — )^cos  r 
2  r0 


(16) 


(17) 

(18) 
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and  the  angle  y'  is  derived  from  Equation  (15)  using  selected  constants  R0, 
kx  and  the  apparent  elevation  tracking  angle,  90. 

At  long  ranges  and  zero  elevation  angle,  Equation  (16)  reduces  to 


AR 


N< 


10c 


(19) 


Using  available  ray-tracing  data  for  the  CRPL  exponential  model  atmosphere1 
for  R  greater  than  2600  km,  and  with 


Ns  =  344.  5 
c  =  0.  157 


then  the  range  error  is  listed  to  be  0.  1164  km. 

Substituting  in  Equation  (19)  gives 
k0  =  0.0189 

Also  for  9o>200  mr  it  was  desired  that 
sin  y  sin  90. 

Again  comparing  with  ray-tracing  data  the  constant  kx  was  selected  to  be 
200;  therefore, 

sin  V  =  sin  0O  +  0.  0189  exp  [-200  0O2]  (20) 

Figure  2  shows  a  comparison  of  calculations  using  Equations  (lb),  (17),  (18), 
and  (20)  with  available  ray-tracing  data 

A  range  of  600  km  was  selected  for  the  comparison  in  order  not  to 
eliminate  the  exponential  terms  in  Equation  (16). 
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THE  EFFECT  OF  REDUCED  RANGES 


As  the  range  decreases  the  exponential  ratio  in  brackets,  Equations 
(8)  and  (9),  approaches  unity  and  causes  the  range  error  and  path  bending 
magnitudes  to  decrease,  which  is  meaningful  However,  at  very  small 
elevation  angles,  the  quantity  g  becomes  small  as  shown  by  Equation  (11). 

Now  as  the  range,  R,  decreases  the  condition  on  (k  +  g),  Equation  (12),  be¬ 
comes  invalidated  This  combination  of  small  angles  and  reduced  ranges 
again  sets  a  condition  which  reduces  the  accuracy  of  the  calculations. 

In  the  same  way  as  sin  was  modified  to  maintain  accuracy  at  small 
angles  a  further  modification  was  used  to  constrain  the  range  dependent  effect, 
again  at  small  angles 

A  convenient  and  useful  form  was  to  replace  V1  by  a  new  angle  where 

sin  y  =  sin  0O  +  [  k0  +  ks  exp  (~k3  R2)  ]  exp  (-kx  902).  (2l) 

Sin  y  then  approaches  the  original  form,  Equation  (15),  as  the  range 
becomes  large  and  approaches  sin  Q0  under  all  conditions  as  the  elevation 
angle  becomes  large.  At  some  minimum  elevation  angle,  for  example,  10  mr, 
sin  y  is  set  equal  to  a  function  of  k2,  k3,  and  range.  Using  two  extreme 
range  situations,  such  as  R  equals  50  km  and  600  km,  tabulated  ray-tracing 
data  will  provide  values  for  the  range  errors  from  which  kg  and  k3may  be 
determined. 

The  simplest  approach  is  to  assume  that  exp  (-k3R2)  is  unity  at 
ranges  under  100  km  and  ks  is  then  directly  determined  to  provide  agreement 
between  the  ray-tracing  data  and  the  new  form  of  Equation  (16),  where  sin  y 
is  now  defined  by  Equation  (21)  The  further  requirement,  having  selected 
kg,  is  that  the  factor 

k2  exp  (-k3R2)  <<k0  (22) 

for  R0  >  600  km . 

This  latter  requirement  does  not  impose  any  difficulty  on  a  suitable  selection 
for  k3  once  k2  is  determined 

Although  the  modification  to  the  initial  equations  has  centered  around 
discussion  of  the  range  errors,  the  bending  equation  is  affected  by  all  the 
same  considerations,  therefore,  the  final  form  of  Equations  (8)  and  (9)  be¬ 
comes  : 
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(23) 


T~  -  — b~® -  [1  -  - g  eXP  ^  1 

10a  tan  y  (k  +  g)  exp  [(k  +  g)  ] 

sin  y -  sin  0O  +  [k0  +  k5  exp  (~k3  R8)]  exp(-ki  0O8) 


(24) 


(25) 


(26) 


(27) 


and  the  angle  yis  derived  from  Equation  (25). 

Figures  3  and  4  show  a  comparison  of  calculated  values  of  range  error, 
Ar,  and  path  bending,  T,  respectively,  with  the  CRPL  tabulated  ray-tracing 
data  The  constants  in  the  sin  y  function  were  selected  to  provide  reasonable 
agreement  over  a  wide  range  of  possible  radar  tracking  conditions,,  Certainly 
in  practice  for  a  particular  region  of  angles  and  ranges  of  interest  the 
selection  of  the  four  constants  could  be  optimized 


A  least  squares  smoothing  program  will  be  available  shortly  to  allow 


a  rapid  calculation  of  all  constants  to  give  the  best  overall  agreement  with 
tabulated  ray-tracing  data 
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5. 


THE  SENSITIVITY  OF  CALCULATED  MAGNITUDES  WITH  RESPECT 
TO  CHANGES  IN  THE  EXPONENTIAL,  MODEL  PROFILE 


The  previous  selection  of  constants  was  related  to  a  CRPL  model  with 
a  surface  value  of  344  5  N  units  As  the  surface  value  changes^  the  decay 
constant,  c,  is  correspondingly  changed  in  accordance  with  the  equation  de¬ 
veloped  for  this  CRPL  model 

A  severe  perturbation  was  introduced  in  the  calculations  by  choosing 
a  station  value  of  Ns  equal  to  200  N  units  and  wi*h  a  corresponding  decay 
constant  of  0  118  Figures  5  and  6  show  calculations  of  the  range  error,  AR, 
and  the  overall  path  bending,  T,  again  compaied  aga  nst  CRPL  ray-tracing 
data  for  the  above  particular  profile  and  for  an  mermed'aie  and  long  range 
case. 


The  difference  between  the  seis  of  curves  has  increased  most  sig¬ 
nificantly  under  two  degrees  (35  mr)  elevation  angle,  The  original  refractivity 
profile,  with  a  station  value  of  344  5  represents  very  moist  conditions  where¬ 
as  the  station  condition  of  200  represents  dry  cond  tions  well  above  sea  level 
This  station  value  could  represent  trackrng  from  an  airborne  platform. 

The  significant  conclusion  is  that  the  equations  can  be  used  to  repre¬ 
sent  a  wide  range  of  possible  tracking  situations  and  environments  In 
practice  one  would  select  the  constants  to  provide  good  agreement  lor  the 
average  local  environmental,  conditions  and  the  equations  could  then  effectively 
absorb  variations  above  the  average  without  one  hav  ng  to  regenerate  a  new 
set  of  constants 

The  CRPL  model  was  used  in  the  prevjous  comparisons  because  tb» 
ray-tracing  data  were  available  in  a  convenient  form  Of  course,  these 
equations  can  be  adapted  to  any  other  exponential  model  providing  a  set  of 
accurate  ray-tracing  data  are  available  to  determine  the  constants 
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A  METHOD  TO  DETERMINE  THE  ELEVATION  ANGLE  ERROR 


In  general,  the  ray-path  bending  calculations  have  little  application 
because  only  at  extremely  long  ranges,  in  fact,  when  the  target  is  outside 
the  real  atmosphere,  does  the  overall  bending,  T,  approach  the  value  for  the 
elevation  angle  error,  *  The  elevation  angle  error  is  very  important  in¬ 
formation  in  most  tracking  situations;  particularly  where  position  informa¬ 
tion  is  required  with  only  a  knowledge  of  apparent  range  and  elevation  angle 

From  Figure  7  the  angle,  91?  is  given  by 

ex  sT+(T  e0)  (28) 


Also,  the  slope  of  the  ray  at  any  point,  P  (x,  y)  is 

dy  7T 

^=tan(r-0J 


,  dx 

and  sin  0X  =  — 


Substituting  (28)  and  (30)  into  (29)  gives 
R 

y  =  J  sin  (90  -  T)  dR 


(29) 

(30) 


(31) 


As  mentioned  previously,  the  geometric  difference  in  length  between  the 
measured  ray- path  range  and  the  slant  range  is  small  From  Equation  (31) 
the  true  elevation  angle,  80l  can  then  be  approximated  by 

o  1  ? 

sin  80~*  —  J  sin  (0g  -  T)  dR  (32) 

o 

Therefore,  the  elevation  angle  error,  €,  is 

1  R 

Ci*  e0  -  sin  1  [—  J  sin  (0O  -  T)  dR]  (33) 

o 
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dy 


FIGURE  7  RAY  GEOMETRY  UvSED  TO  CALCULATE  THE  ELEVATION 
ANGLE  ERROR 
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Replacing  T  from  Equation  (24)  leads  to  a  solution  since  the  integral 
is  then  a  function  of  R  However,  the  term  tan  y  in  the  denominator  of  (24) 
complicates  the  integral  and  at  this  time  a  solution  has  not  been  found.  The 
second  approach  was  to  remove  tan  7  from  the  integral  by  considering  its 
mean  value  over  R  This  approach  led  to  significant  errors  at  small  eleva¬ 
tion  angles . 

The  most  satisfactory  solution  was  to  evaluate  the  integral  directly 
by  the  summation  of  terms,  where 

R  n 

J  sin  (0O  -  t)  dR  ~  £  sin  (0O  -  |t|.)  -  (6^)  (34) 

o 

where  Tx  is  the  value  in  the  1 ^  interval 
and  6  Rj  is  the  i ^  interval. 

Figures  8  and  9  show  a  comparison  of  results  for  two  ranges  of  600 
and  200  km,  respectively,  and  with  the  integral  approximated  by  ten  terms 
in  the  summation  The  difference  between  results  is  most  noticeable  at  the 
short  ranges  However,  referring  to  the  bending  calculations,  Figure  4,  it 
is  seen  that  the  T  values  are  not  in  exact  agreement  with  the  ray- tracing 
data;  therefore,  one  would  not  expect  the  €  values  to  be  any  better. 

Since  the  values  for  |t  |j  are  not  linearly  dependent  upon  range  con¬ 
sideration  should  be  given  to  a  weighting  factor  which  would  increase  the 
number  of  terms  in  the  sum  at  shorter  ranges.  This  approach  is  being  con¬ 
sidered  in  the  optimization  analysis  and  results  will  be  presented  at  a  later 
date . 


At  long  ranges,  T  becomes  essentially  independent  of  range  and  the 
integration  of  Equation  (33)  gives 


C  =  0O  -  sin*"1  [sin  (0O  -  It  |)J 


(35) 


In  this  case,  Equation  (35)  gives 

€=  M  (36) 


As  ray  tracing  data  showsi?  this  approximation  is  only  useful  where 
tracking  ranges  in  excess  of  3000  km  are  involved 
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AN  APPROXIMATE  SOLUTION  FOR  THE  DOPPLER  VELOCITY 
ERROR  ANGLE 


Referring  to  Figure  10,  the  doppler  velocity  ei  ror  in  the  radial 
direction  attributed  to  refraction  can  be  shown  to  be 

A  V  ~  V  sin  il>  (  6)  (37) 

where  V  is  the  actual  target  velocity 

t /)  is  the  angle  between  the  velocity  vector  and  the  direct  slant  path,  R0, 
and  6  is  the  local  angle  between  the  ray  direction  and  the  slant  path 

From  the  measurements  of  the  range,  apparent  elevation  angle,  and 
therefore,  the  calculated  elevation  angle  error  (Equation  (33)),  it  is  possible 
to  determine  the  apparent  target  velocity  and  its  direction  of  motion  Then 
the  calculation  of  the  angle,  6,  will  provide  a  first  order  correction  to  the 
velocity  measurement  error  produced  by  the  ray  path  bending. 


From  Equation  (32),  sin  $0is 

sin  ^0-  \ 

R 

J  sin  (0O  -  T)  dR 

.  38) 

o 

Expanding  sin  (90  -  T)  and  assuming  T  to  be  a  small  quantity  which  is 
generally  true,  then 


sin  (9o*  T)  ~  sin  0O  -  T  cos  Q0 


(39) 


Substituting  (39)  in  (38)  and  substituting  for  the  magnitude  of  T  f  r  om  (24)  gives 

R 

Q  .  „  Ns  cos  0O  (.  dR 

s,„  S,~  sm  8,  -  — 

o 


Ns  cos  0n 
+  - 


R 


J  s 


g  exp  [&fl 


10°R  ^  tan  V ( g  +  k)  exp  [\g  +  k)2J 


dR 


(40) 


For  long  ranges  v  is  essentially  independent  of  range,  according  to 
the  method  of  selecting  the  constants  given  by  Equation  (22)  Then  for  small 
elevation  angles,  Equation  (40)  can  be  written 


Ns  cos  0o  N s  cos  9 n  n  g  exp  [g2]  dR _ 

10  tan  v  10  R  tan  y  '  (k  +  g)  exp  [tk  +  g)  j 


(41) 
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Vr  =  V  COS  (*  +  8) 


FIGURE  10  DOPPLER  ERROR  ANGLE  GEOMETRY 
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From  the  definition  of  g  and  k,  Equations  (26)  and  (27),  the  last  term  in 
Equation  (41)  can  be  written 


?  B  exP  1  dR  _  ?  exp  [-  (a  R  +  2  a  g  R)]  dR 

0  (k  +  g)  exp  [(k  +  g)2]  J  (a  r  +  g) 


i  ,  C  v- 

where  a  =  (- -  )  cos  V 

2  r o 

Completing  the  square  in  the  exponent  and  setting 
x  =  a  R  +  g 

then  Equation  (42)  becomes 

f  g  exp  [g2]  dR  _  g  exp  [g2]  (aR+g>  exp  [-x2]  dx 

“  (k  +  g)  exp  f(k  +  g)  ]  a  -  x 

S 

r  2-1  taR+g)2 

_  g  exp  [g  J  n  exp  t-t) 

2a  i  t 

g 

where  x2  =  t 


(43) 


(44) 

«  45) 


The  integral  of  Equation  (45)  can  be  evaluated  by  expanding  and  integrating 
to  give 


oo 


,ap+g)  exp  (-t)  dt.  =  ln  [Jl±j.]+  s  Liyi  <k  +  g) 

2  t  g  n=i  n.  n! 


2n 


-  r 


00  (-Dngsn 


n=i  nu  n ! 
2 


(46) 

(47) 


Ei  [g2J  -  Et  [(k  +  g)2J 

where  E^(z)  is  the  exponential  integral  defined  by 
00 

*lW  ■  I  *> 

z 

Then,  from  Equation  (41),  ^Qcan  be  written 

- Tfcct°£nv° + [t‘]  ie. [gs]-E. f<k+s>EJj  <«> 


10  2  k  tan  y 


where  k  =  a  R 
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From  Figure  1, 


0o-  9  o  “  € 

=  90  -  T  +  6 


(49) 


A  comparison  of  Equations  (24),  (48),  and  (49)  suggests  that  at  long  ranges 
then  Tcan  be  identified  with  the  second  term  in  Equation  (48)  at  small  eleva¬ 
tion  angles  Also,  from  (48)  and  (49)  the  doppler  velocity  error  angle,  6, 
can  be  identified  with  the  last  term  in  (48),  and 


5~  Ns  cos  6n  g  exp  [gs] 
10s  2k  tan  V 


[Ei  [g2]  -  EA  [(k  +  g)2]] 


(50) 


Figure  11  shows  a  comparison  of  the  calculated  value  of  6  with  ray¬ 
tracing  data  A  range  of  600  km  was  chosen  to  ensure  that  V  would  be 
independent  of  range. 

The  agreement  is  reasonable  except  at  very  small  angles  where,  of 
course,  the  previous  calculations  for  €  also  deviated  from  the  ray-tracing 
data. 


The  validity  of  the  solution  is  again  dependent  upon  having  a  range 
sufficiently  large  such  that  V  is  range  independent  and  that  the  elevation 
angles  are  sufficiently  small  to  justify  the  expansion  in  Equations  (39)  and 
(41).  Since  the  angle,  6,  does  not  become  significant  in  practice,  except  at 
lower  elevation  angle  tracking,  the  latter  condition  does  not  impose  any 
practical  difficulty. 
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8. 


SOME  FINAL  COMMENTS 


It  is  evident  in  the  development  of  the  elevation  angle  error  equation 
that  the  calculation  of  e  will  be  no  better  than  the  ray-path  bending  values  for 
T  In  turn,  the  €  values  directly  influence  the  values  obtained  for  6,  the 
doppler  velocity  error  angle.  In  practice,  therefore,  if  one  is  satisfied  with 
using  an  exponential  model  to  define  the  propagation  conditions,  every  effort 
should  be  made  to  select  the  k  constants  to  give  good  calculations  for  T  over 
the  range  of  interest. 

Referring  to  Equations  (23)  and  (24),  it  is  seen  that  the  range  error 
and  the  path  bending  are  related  by  the  equation 

T  =  -  c  cos  r  (AR).  (51) 

For  real-time  tracking  operations  this  simple  relationship  provides  immediate 
values  of  path  bending  as  the  range  errors  are  being  calculated  In  practice, 
this  ability  is  significant. 

The  basic  development  of  the  equations  has  suggested  other  approaches 
and  refinements  which  were  not  pursued  in  this  paper.  The  attempt  herein 
was  to  derive  the  simplest  types  of  equations  for  use  in  real-time  tracking 
applications  at  the  expense  of  rigor  and  maximum  possible  accuracy.  The 
development  of  a  least  squares  smoothing  program  will  provide  a  method  to 
determine  the  constants  k0  through  k3  to  give  the  best  comparison  with  any 
particular  set  of  ray-tracing  data.  Two  ray-tracing  programs  presently 
available  can  be  used  to  generate  tabulated  data  for  any  exponential  model 
atmosphere . 
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